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Abstract. We give a short review of known exact inequalities that can be 
interpreted as ‘energy—time’ and ‘frequency-time’ uncertainty relations. In 
particular we discuss a precise form of signals minimizing the physical frequency¬ 
time uncertainty product. Also, we calculate the ‘stationarity time’ for mixed 
Gaussian states of a quantum harmonic oscillator, showing explicitly that pure 
quantum states are ‘more fragile’ than mixed ones with the same value of 
the energy dispersion. The problems of quantum evolution speed limits, time 
operators and measurements of energy and time are briefly discussed, too. 


This is our present to Margarita Alexandrovna and Vladimir Ivanovich Man’ko 
on occasion of their 75th birthdays. 


1. Introduction 

The energy-time uncertainty relation (ETUR) 

AEAT > h (1) 

is one of the most famous and at the same time most controversial formulas of quantum 
theory. It was introduced by Heisenberg [I] together with his coordinate-momentum 
uncertainty relation [|] 

Ax A p > h. (2) 

The importance of both relations for the interpretation of quantum mechanics was 
emphasized by Bohr [2], However, the further destiny of relations (JTJ) and ([2]) turned 
out quite different. A strict formulation of relation ([2]) was found almost immediately 
by Kennard 3 and Weyl [2] in the form of inequality 

AxAp > H/2, (3) 

where Ax and A p are well-defined quantities, namely mean-square deviations 
(dispersions). A liitle later, more general inequalities were derived by Robertson [5j[6] 
and Schrodinger |7 . Actually, the meaning of inequality © is rather different from 
the original thoughts of Heisenberg and Bohr related to measurements of canonical 
variables (see e.g. [8] 111). Nonetheless ([3]) is a well established consequence of the 
quantum mechanical formalism. 

$ Heisenberg and Bohr used the symbols = or ~ instead of >. 
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On the contrary, the physical and mathematical meanings of inequality ([I]) 
appeared to be much less clear than that of § (and less clear than Heisenberg, 
Bohr and other creators of quantum mechanics thought initially). The main reason 
is that, in fact, there are several quite different physical problems where relations like 
([I]) can arise, and in each concrete case the meaning of the quantities standing on the 
left-hand side proves to be different. This was clearly demonstrated for the first time 
by Mandelstam [§] and Tamm 12 and by Fock and Krylov 13 , and many authors 


arrived at the same conclusions later 14-17 


The aim of this paper is to give a concise review of some results, where the 
strict inequality sign (>) replaces the symbol > of approximate (by an order of 
magnitude) inequality in relation ([Tj) . Such a mini-review seems to suit well for 
the inclusion to the special issue celebrating 75+75 years of Margarita and Vladimir 
Man’ko, because they have been always interested in the related subjects: see, 
e.g., references 18 -20 where different aspects of time-energy and time-frequency 
uncertainty relations were discussed (together with many other topics). Our initial 


18 19 


goal was to extend or generalize some results obtained many years ago in 
These generalizations are given in sections m However, it seems reasonable to 
describe the most important achievements of many other authors as well, providing 
the relevant references. Although many reviews on the ETUR were published during 


decades 21-25 , we hope nonetheless that this mini-review with the updated reference 


list will be also useful for many readers as the complementary source of information. 


2. Mandelstam Tamm inequalities and their application to the problem 
of decay 

The first rigorous formulations of relation |T]) was given by Mandelstam and Tamm 
12 . They used the uncertainty relation for two arbitrary Hermitian operators A and 
B (derived for the first time by Robertson hj) 


AAAB > - 
~ 2 


dAB}) 


(4) 


Then, choosing B = H (where H is the system Hamiltonian) and remembering that 
operator (i/h)[H, A] is the operator of the rate of change of the quantity A, i.e., 
(i/H)[H,A\ = dA/dt (provided operator A does not depend on time explicitly), one 
can rewrite Q as 

AEAt A > h/ 2, (5) 

where AE = AH and 

At a = AA/\(dA/dt)\. (6) 

The meaning of inequality ([5]) is that it yields an estimate of the time interval required 
for a significant change in the average value of observable A: by an amount of the 
order of the mean squared variation. It results in the statement that ‘a dynamical 
quantity cannot change, remaining always dispersionless’ 12 . 

Relation ([Hj) may seem a little artificial at first glance, since operator A may be 
quite arbitrary. However, there exists at least one important specific choice of this 
operator. Namely, let us consider, following 12 , the projector L on some initial 


§ We use the name Mandelstam as it appears in western literature; in Russian literature it is written 
as Mandelshtam. 
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quantum state |^>(0)): L = |-i/(0))(V’(0)|- Since L is the projection operator, L = L 2 
and 


Then inequality ([5]) assumes the form 

< A E^L-(Ly. (7) 

Integrating ([ 7 ]) with account of the initial condition L(0) = 1, Mandelstam and Tamm 
obtained relations 


ft dL ft 
2 ~dt ~ 2 


where 


A Et/h > 7t/ 2 — arcsin \JQ(t), 0 < t < 00 , 

(8) 

Q(t) > cos 2 (A Et/h), 0 < t < ttH/2AE , 

(9) 


(10) 


One may interpret function Q(t) as the probability to remain in the initial state 
IV'(O)). In this case it seems natural to define the half-decay period Tj/ 2 by means of 
the relation Q (Tj/ 2 ) = 1/2. Then (j9| results in the inequality 

T 1/2 AE > tt/j/ 4 « 0.785ft (11) 

with well defined quantities T]/ 2 and A E. 

An immediate important consequence of inequality <§ is the impossibility of 
strictly exponential decay 


Q{t) = exp (—t/r) 


( 12 ) 

Indeed, 


for realistic physical systems with a finite energy dispersion A E [T2|jT3| 26,27 
it follows from (|9| that the law (12) can be realized only approximately and for 
sufficiently big values of time satisfying the inequality 

t > r In [1 + ft 2 /(2rAfti) 2 ] . (13) 

Generalizations of relations ([8]) and to the case of time-dependent 
Hamiltonians were obtained in 28 . Noticing that inequality (|5 1 becomes meaningless 


in the case of A E = 00 , Mandelstam and Tamm wrote in 12 
desirable to find a more general relation of the same type as (15 
sections we demonstrate how this can be done and what concrete generalizations were 
proposed. 


that ‘it would be 
In the following 


3. Decay laws and spectral distributions 


The relations between the decay law Q{t) and the energy spectrum of the system were 
established by Krylov and Fock [l3]. Suppose that we know the decomposition of 
vector |V’(0)) over the energy eigenstates: 

1^(0)} = J a(E)\E)dE, (E'\E) = S(E — E'). (14) 


Then 


IV’O)} = 


a(E) exp (— iEt/H ) | E) dE 


(15) 
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and the nondecay amplitude y(t) = (^(0)|^(f)) (called ‘integrity amplitude’ in [26] ) 
can be expressed as the Fourier transform of the positive energy distribution function 
P(E) = |a(£)| 2 : 


x(t) = (V’(O)I = J P{E) exp (—iEt/h) dE. 


Its consequence is the identity 

J P{E)dE = 1 


(16) 


(17) 


The probability of finding the system in the initial state at time t equals 
Q(t) = \x(t)\ 2 - The energy variance can be calculated as 


with 


(A E) 2 = J (E - E) 2 P{E)dE 
E = I EP(E) dE. 


(18) 


(19) 


Fock and Krylov 13 proved that the necessary and sufficient condition of decay (the 
property Q(t) —> 0 for t —t oo) is the continuity of the integral energy distribution 
function P(E) = f E P{e)de. This means, in particular, that the energy spectrum must 
be continuous, in order that function P(E) would not contain terms like S(E — Eq) 
corresponding to discrete energy levels. 


It is known that the exponential decay law (12) corresponds to the Lorentzian 
energy distribution 

F/2tt 


P(E) = 


( 20 ) 


(E — Eq) + T 2 /4 

The peak width T, defined as the length of the energy interval where P(E) > P{ 0)/2, 
is related to the lifetime r by the equality 

rT = H. (21) 


However, distribution (201 is an idealization (although very good one in many practical 


cases), because it results in the relation A E = oo. Another drawback of distribution 
20l is that it implies that the energy spectrum stretches from —oo to oo (only under 


this assumption integral (16) yields the exponential function of time). But the energy 
of real physical systems is limited from below, and this fact leads to violations of the 
exponential decay law for t —> oo, when Q(t ) ~ t~ 13 with some constant ft depending 
on the concrete form of the energy spectrum [29] (see [27 30-34] for later discussions 


and reviews). Moreover, an oscillatory decay is also possible 


35 


A simple example of the ‘decay’ that never has the exponential form was given 
a freely expanding Gaussian wave packet in one dimension has 


by Bhattacharyya 36 


the energy distribution function 

_*j^ j 2 

P(E) = ( 'V2irEAE ) exp (-E/V2AE^J 

and the nondecay probability 

Q(t) = (1 + 2(A E) 2 t 2 /h 2 )~ 1/2 . 

In this case T 1/2 AE = 3/2. 


( 22 ) 


(23) 


















Energy-time and frequency-time uncertainty relations: exact inequalities 


5 


The remarkable formula (16) has many interesting consequences. One of them was 
obtained by Luo [37 . We reproduce it in a slightly modified form. Let us consider 
function A(y;t) = (Re [exp(iy)x(f)]) 2 , where y is an auxiliary variable. Obviously 
Q(t) = |x(t)| 2 is the maximal value of A(y\t) for the fixed value of t. On the other 
hand formula (16) yields (we put h = 1 here) A{y\t) = [/ cos(y — Et)P(E)dE] 2 . 
Using the Schwartz inequality 

r- 2 


\f{x)\ 2 dx / \g(x)\ 2 dx > 


Re [f(x)g*(x)\dx 


(24) 


with x = E, f(E) = cos (y — Et)-^P(E) and g(E) = y/P(E) we get, taking into 
account ©> 

A(y, t)< j cos 2 {y - Et)P(E)dE = ^ j [1 + cos(2 y - 2 Et)} P{E)dE 

\xm\- 


= \ (1 + Re [exp(2iy)x(2t)]) < ^ [1 


Finally we have 


Q(t) < - [l + y/Q(2t) 


(25) 


Some interesting consequencies of this inequality were discussed in 37 . The 


importance of inequality (25) consists in the fact that it forbids many decay laws 


that one could invent ‘from a head’. For example, it clearly forbids the instantaneous 
decay, such that Q{t) = 1 for t < f* and Q(t) = 0 for t > f*. However, it is not 
the strongest possible inequality, because it does not forbid the exponential decay law 
(121, which is forbidden by § lor short times. On the other hand, looking at (|9| one 
could suppose that the none-decay probability could behave as Q{t) ~ 1 — at° with 
b > 2 for t —> 0. But inequality (25) can be fulfilled in this limit for b < 2 only. So, 


only a parabolic time dependence of Q{t) is permitted for short times. This is clear 
from the well known short-time Taylor expansion 


Q(t) = 


(V’olexp - 


iHt 


IV’o) 


1 - 


it(H) 

h 


t 2 (H 2 ) 
2 h 2 


= 1 - 


h 2 


((iL 2 ) - ( H ) 2 ) (26) 


4. Wigner’s approach 


Wigner 38 studied the function \ u {t) = (u\ip(t)), i.e., the probability amplitude of 


finding the quantum system described by the vector | ip(t)) in some arbitrary fixed 
state |u). The Fourier transform of function Xu{t) was defined as (let us put here 
H= 1 ) 


y(E) = (2t r) 1/2 J x u {t) exp (iEt)dt. 


(27) 


If the energy spectrum (i.e., the spectrum of the Hamiltonian operator) is absolutely 
continuous and bounded from below, then, taking E m i n = 0, one can expand vector 
|u) over eigenstates | E) of the Hamiltonian operator as follows, 


u) = 


b{E)\E)dE. 


(28) 


Taking into account equation (15) we have 


r)(E) = (2Tr)^ 2 b*{E)a{E). 


(29) 
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Wigner defined the characteristic time tw of the system’s stay in the state |u) and 
the energy dispersion £\y as follows (assuming that the initial time is zero): 


•w 


~2 

£ W 


ir t 2 \xu(t)\ 2 dt 
ir\xu{t)\ 2 dt 1 
^(E-Eo^E^-dE 
J 0 °°\r,(E)\ 2 dE 


(30) 


(31) 


E 0 being an arbitrary parameter. If 77 ( 0 ) 7 ^ 0, then function Xu(t) behaves as 1/t for 
t —> 00 and t = 00. Therefore Wigner considered the states satisfying the restriction 
77 ( 0 ) = 0. He obtained the formula 


r 

_2 _ Jo 
T w ~ 


\dr]/dE\ 2 dE 
f 0 °° \ V (E)\ 2 dE ' 


(32) 


Using the Schwartz inequality (24) he arrived at the inequality (recovering the Planck 
constant) 

£w T w > h/ 2. (33) 


The equality in (33) cannot be achieved due to the boundedness of the energy spectrum 
and the restriction 77 ( 0 ) = 0. (The equality £w T w = 7i/2 holds for the Gaussian 
distributions; but these distributions obviously do not satisfy the imposed restrictions.) 
The minimal possible value of the product £w T w depends on the parameter Eq. For 
E 0 = 0 this minimal value equals [38] ( £w T w)min° = 37i/2. It is achieved for the 
function t]q{E) = E exp[—3-E 2 /(4e 2 )]. 

A comparison of the M T and Wigner UR was performed in 


16 for free one¬ 


dimensional wave packets. It was shown that both methods lead to almost identical 
results. For the recent analysis of the Mandelstam-Tamm UR one can consult 39 . 


5. Different decay times of unstable systems 


Besides the half-decay time T)/ 2 and the Wigner time tw (30), many other definitions 
of the decay time are possible. Fleming 26 suggested to use the quantity 


To = 


Q(t) dt 


(34) 


where the non-decay probability Q(t) was defined by equation (10). Definition (34) 
gives the lifetime tq = r for the exponential function (12). Taking into account 
inequality Q and writing z = ttH/(2AE) we obtain the following lower bound for 
TO 


26 


r 0 > f Q(t) dt> f cos 2 (A Et/h) dt = 
Jo Jo 


irh 

4A E' 


(35) 


A weaker inequality tqAE > h/2 was obtained in [40] 


The most strong inequality with an achievable lower bound was given in study 41 


First we note that the consequence of equations (10), (16) and (341 is the formula 


tq = 7 rh 


dE[P(E)} 2 . 


(36) 
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Therefore one has to minimize the functional 


/ oo r roc 

dE[P(E )] 2 / dEP(E)E 

-OO \_J —oo 


1/2 


where the origin of the energy scale is chosen in such a way that E = 0. The value 
of functional I(P) is not changed if instead of function P(E) we use the function 
P\(E) = XP(XE) with an arbitrary (positive) parameter A. Choosing this parameter 
in such a way that (E 2 ) = 1, we arrive at the problem of minimizing the functional 


m = 


dE[P(E)f 


(37) 


under the following auxiliary conditions (here we use dimensionless variables and 
ft = 1): 


/ OO fOO 

dEP(E) = / dEP(E)E 2 

-oo J — OO 


= 1, 


(38) 


(39) 


/ dEP(E)E = 0, P(E) > 0. 

J —OO 

It was shown in [41] that the extremal function is the truncated parabola 

/ (n/45/20) (1 - E 2 / 5 ) , \E\<y/5 

P ” (E) -\ 0, |E|>^5 ■ 

To prove this result, let us consider some function R(E) satisfying conditions (38) and 


(40) 


([39]). If f(E) = R(E) - P 0 (E), then 

I(R) = I(P 0 + f) = I(P 0 ) + 1(f) + 2 


dEP 0 (E)f(E). 


(41) 


It is obvious from (37) that 1(f) > 0; moreover, 1(f) >0 if the function f(E) can 


turn into zero only for a set of zero measure. Let g(E) be the non-truncated parabola 


(40). Since R(E) and Pq(E) satisfy conditions (38) and 

/»oo 

/ dEf(E)E n = 0, n = 0,l,2. 


Consequently, 


rV 5 


0 = / dEf(E)g(E) = / dEf(E)g(E) 

J—OO J—y /5 


'I E\>V5 


dEf(E)g(E). 


(42) 


(43) 


The last integral in (43) is nonpositive, since g(E) < 0 and f(E) = R(E) > 0 for 
\E\ > y/5. Consequently, the first integral on the right-hand side of (43) [which is 
identical to the last integral on the right-hand side of i s no nnega tive . Therefore 
I(R) > I(Pq) for any function R(E) satisfying conditions (38) and (39). Thus we 
obtain the inequality 


. 3tt _ 
T ° AE -^ k ' 


0.843ft. 


(44) 


The equality sign in (44) is achieved for the distribution P*(E) = (AE)~ X Pq(E/AE) 
(here we return to dimension variables), where Po(E) is given by (40) and the energy 
dispersion A E can be arbitrary positive number (remember that the energy scale is 
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shifted in such a way that (E) = 0). The corresponding ‘nondecay probability’ can be 
calculated with the aid of formulas (101 and (16). It has the form 


Q*(t) = 9 (sin 2 — .zcosz) 2 /z 6 , (45) 

where z(t) = \/5tAE/h = 37rt/(5ro). 

Two other examples of energy distributions possessing products tqAE close to 


the minimal possible value (|44j) were given in 141 . The first of them is the Gaussian 
distribution 


Q(t) — exp 7r (f/2ro) 


= exp 


- ( tAE/ny 


P(E) = (A eV2tt 


-l 


exp 


~\ ( E/AEf 


t 0 AE = v^F/2 « 0.886. 


(46) 


(47) 

(48) 


One can check that the nondecay probability (46) satisfies the Luo inequality (25). 


The second example is the stepwise energy distribution 

\/3/(6A E), \E\ < VZAE 

0, \E\ > \/3A E 


P{E) = 


Q(t) = 


sin (V3tAE/ti) 
V3tAE/h 


(49) 


(50) 


t 0 AE = —= « 0.907. 
2y 3 


The time-dependences (45), (46) and (50) seem to be very far from those which can 


occur in real decaying physical systems. However, they can be interesting from the 
point of view of the problem of maximal speed of quantum evolution: see section [8] 
Mathematical aspects of the time decay problem (time asymmetry) were discussed 
in [421446 . 


5.1. Modifying definitions of the decay times and energy spread 

The energy distributions in realistic decaying systems are close to the Lorentz 
distribution (20). In these cases the decay time is determined not by the energy 
dispersion ©7 but by the energy level width T. Changing the form of the distribution 
function P(E) at its ‘tail’ (for | E — Eq\ T) we can change the variance A E 
significantly, but the decay time will remain practically unchanged. This fact was 
emphasized long ago by Fock and Krylov [13 . Therefore, to fill inequality 0 with 
a physical content in the decay problems, we must find a more reasonable definition 
of energy ‘uncertainty’ A E (as was questioned by Mandelstam and Tamm), linking 


it not to the energy variance (18), but to some other quantity, in such a way that it 


would be close to the energy level width for distributions similar to the Lorentz one. 

Several possible definitions of this kind were proposed in [19] . They were based 
on some results of study 16 , where the concept of an ‘equivalent width’ of a function 
was introduced. Precisely, the equivalent width W{<p) of function tp{x) was defined as 

/ OO 

y{x)dx/w{Q), (51) 

-OO 
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provided the integral exists and <^(0) 7 ^ 0. Further calculations are based on the 
following simple observation: if two functions f{x) and f(y) are related by the Fourier 
transformation, i.e. 


then 


M = / e~^f{x)dx, 


W(f)W(f) = 2tt. 


f(x) = / e lxy f(y)dy/{2n), 


(52) 


Now let us look at the following consequence of equation (16 1 (in dimensionless 
variables with h = 1 ): 

/ oo 1 roo 

e~ iEt [P(E )] 2 dE = — J- x ( t ') X * (f - t) dt 1 , 

where P(E) is the energy distribution function and x(t) is the non-decay probability 
amplitude. Defining f(E) = [P{E + Eg)] 2 , where Eo is an arbitrary real number, we 
have 


Obviously, 


whereas 


Obviously 


f(t) = exp (iE 0 t) / x(t' + t) X* (*') dt'/(2Tr). 

J — OO 

/ OO 

[P{E)f dE/[P(E 0 )] 2 , 

-OO 

/ / dtdtf exp (iE 0 t)x (' t ' + t) x* {?) 


W(f) = 


mf)\ < 


I-o o lx (t) I 2 ) dt' 

f f dtdr\ X (r + t) x* (t) | _ [/ dr\ X (r) |] ‘ 


fZo Ix(t) \ 2 )dr 


f\x(r) 1 2 )dr ' 


In addition, the identity \W(f)W(f) \ = 2n holds for any value of E 0 as a consequence 
of (52). Remembering that x(—t) = X* (t) for real energy distribution function P{E) 
and returning to dimension variables and the nondecay probability Q(t) = |x(^) 1 2 5 we 
arrive at the inequality 

2 


/ I P(E)} dE 


r vwidt 


> 7 tH, 


(53) 


\P(E 0 )} 2 / 0 °° Q{t) dt 

which holds for an arbitrary value E 0 , in particular for E 0 corresponding to the 
maximum of function P(E). Looking at the left-hand side of (53), it seems reasonable 
to introduce the following definitions of the decay time and energy uncertainty: 


T* = 


r dt 


4 / 0 “ Q{t)dt 


Afi, = 


fjpm 2 dE ' 

max [P(E)] 2 


(54) 


Then r* = r for the exponential decay (12) (the normalization factor 1/4 in the 
definition of r* is chosen just in order to ensure this equality). Therefore the lower 
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bound for the product r* AP* is the same as in © or ( |35[ > , but with different meanings 
of symbols on the left-hand side: 

T*Ai?* > nh/A. (55) 


On the other hand, taking into account (34) and (361, we can define the characteristic 
decay time and energy uncertainty as 
1 


/•OO _ 

/ VQ(t) dt, AP** = [maxP(£’)] _ 

Jo 


(56) 


rewriting (53) in the form 

t**A.E** > nH/2 (57) 

(the coefficient 1/2 in the definition of r** is chosen again in order to ensure the 
equality r*» = To for the exponential decay law). Three sets of decay times and energy 
uncertainties are connected as follows: 

T* = t*„/t 0 , a P* = AE^To/nh. (58) 

Note that inequality (57) can be obtained directly from (52) if one takes / = P(E). 
Then / = x(t) according to (16), so that (57) follows from (52) due to the inequality 
W(f) < W(\f\) (the value W(f) is real and positive in this case due to (52), since 
function / = P(E ) is positive). 


distribution (201. The same is true for relation fl55to due to fl58b . Therefore, inequalities 


Relation (57) becomes an equality for the exponential decay law (12) with energy 


(55) and (57) can be considered as reasonable (and exact) energy-time uncertainty 


relation for decaying systems. 

Taking |it) = |^(0))) in the Wigner relations and ( |33[ ), we can rewrite 

them in terms of functions P{E) and Q{t) as follows (relaxing the restriction of the 
energy boundedness from below): 


, 2 _ / 0 °° t 2 Q(t)dt 
fo°° Q(t)dt 


1 f 00 

= — t 2 Q(t)dt, 
to Jo 


e 2 = 


S-^E ~ (E)) 2 [P(E)] 2 dE 


S-jP{E)] 2 dE 


(E) = 


f^°oo E[P(E)) 2 dE 

r[p(E)) 2 dE 5 


st > h/2. 


(59) 

(60) 
(61) 


For the exponential decay law (12) with the Lorentzian energy distribution (20) we 
have 

f = V2H/T : (E) =E 0 , £ = T/2, ts = h/y/2. 

As one can see, the replacement of the weight factor P{E) by [P(E)] 2 in the 
formula for the ‘effective variance’ s 2 enables us to suppress the slowly decreasing 


‘tail’ of the Lorentz distribution function P(E) (20). Asa result, the ‘effective variance’ 


proves to be of the order of the physically acceptable width of energy level. Moreover, 


for the Gaussian decay law (|46j)-( 47) the inequality (61) becomes strict equality 


(similar to the minimal product AxAp = h/2 for the Gaussian distributions in the 
coordinate and momentum spaces). Indeed, we obtain in this case t = h/(y2AE) 
and £ = AE/y/2. 

Fujiwara [47 introduced one more characteristic decay time 

/ 0 °° tQ(t)dt 


Ait = 


n = 


r Q{t) dt ’ 


(62) 
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where t is defined by equation (591. Then we have for the exponential decay law (12) 
ti = Ai = H/T, eAi t = h/2. (63) 

However, the value H/2 is not the lower limit for the product eAi t (as was mentioned 


(461-(47): 


without proof and details in 47 ). It is clear from the example of Gaussian states 


Ti = H/(AE^/n), AH = (H/AE)^/(n - 2)/(2ir), 
(eAi t) G au SS =h[(n- 2)/(4tt)] 1/2 « 0.301A 


(64) 


Other definitions of the ‘evolution time’ were considered in 48 . Concrete calculations 


of functions like f u (t) and ij(E) (27), as well as the related quantities f, e and ri, for 
different choices of the reference states |u) and solutions to the Schrodinger equation 
ip(t) for a partiele in a uniform electric field (when the spectrum of the Hamiltonian 


is continuous and extends from —oo to +oo) were made in study 49 


Hilgevoord and Uffink |50.51] introduced a measure for the uncertainty in energy 
W a (E) as the size of the shortest interval W such that 


\(E\i/j(0))\ 2 dE = a. 


(65) 


iw 


The time extension Tg was defined as the minimal time it takes for IV’(O)) to evolve to 
a state | if(r)) such that 


|(V>(0)|V-(t))| = /?. 


It appears that 


rg \V a > 2tiarccos 


/3 + l — a 


a 


for (3 < 2a — 1. 


( 66 ) 


(67) 


In particular, one obtains for the Lorentzian distribution (20) the values 
W a = T tan(a7r/2), rg = 4r ln(l//3). 


6. ‘Time frequency’ uncertainty relations 

It appears that the problem of finding the lower limit for the product eA\t is closely 
related to the time-frequency uncertainty relation. Therefore it is worth discussing 
this problem in details. 

Dividing both sides of inequality |l]) by the Planck constant h and using the 
Planck formula E = hu> one arrives at the approximate inequality 

AuAt > 1 (68) 

where the right-hand side can be replaced with the same accuracy by tt or 2n. Actually, 
this inequality, which connects some effective duration of signal At with some effective 
spectral width Aca, was known in optics and radiotechnics long before the birth of 
quantum mechanics. Moreover, it was frequently used for illustrations of the quantum- 
mechanical uncertainty relations at the early years of quantum mechanics. However, 


exact formulations of relation (68) were given after strict formulations of the quantum- 


mechanical UR only. By analogy with the strict Heisenberg-Weyl UR AxAp > h/ 2, 


some people believe that the precise version of (68) is the inequality 
Aw At > 1/2, 


(69) 
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where the quantities Aw and At are defined in the same standard way as in quantum 
mechanics: 


/ OO pOO 

0, t= / tf(t)dt , (70) 

-oo J — OO 

/ OO 

(w — w) 2 |F(w)| 2 dw, (71) 

-oo 

/ oo 

w|F(w)| 2 dw, (72) 

-oo 

l r°° 

fit) = -f= / (73) 

V "7T J —oo 
1 /*°° 

F(w) = — / f(t)e~ %ut dt. (74) 

V "TT J —OO 


It is supposed here that the signal /(f) is normalized according to relations 

/ oo poo 

/ 2 (t)dt= / |F(w)| 2 )dw = l. 

-OO J — oo 


(75) 


In such a case, although inequality (|69j) is correct, it is useless in many practical 
situations. The problem is that, in contrast to the quantum-mechanical case, the 
time-dependent signal /(f) is always a real function. Therefore 

F(u>) = F*{-u>), |F(w)| 2 = |F(—w)| 2 , (76) 

so that w = 0 for any signal. This result contradicts our intuition. For example, for any 
narrow-band signal, function |F(w)| 2 has two narrow peaks, centered symmetrically 
with respect to point w = 0. But in such a case the quantity Aw, defined by equation 
(721, characterizes not the spectral width of the signal, but its carrying frequency. 


Therefore, in order to have a meaningful UR for real signals, it seems much 
more natural to take into account only positive frequencies and redefine the average 
frequency and spectral width as follows: 


pOO 

(Aw + ) 2 = 2 / (w - w + ) 2 |F(w)| 2 dw, 
Jo 


w + = 2 / w|F(w)| 2 dw, 

Jo 


(77) 


(78) 


where the factor 2 appears due to the normalization. Thus we arrive at the problem of 
finding the minimal possible value of the product Af Aw+. This problem was analyzed 
long ago (in 1934) by Mayer and Leontovich 52 . However, this study was practically 


unknown for a long time. More known was the paper by Gabor 53 , where different 


definitions of the frequency and duration of signals were considered. Gabor introduced 
the concept of ‘analytical signal’: a complex function /+(t) whose Fourier transform 
does not contain negative frequencies, 


pOO 

f+(t) = / F(w)e^dw/-v/7T, 

Jo 

F+M = f f+^e~ iult dt/V2n = | ^ (w) ’ " ^ ° 


(79) 


(80) 
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Putting the function F + (u>) in the definitions (721, we have the equalities uj = uJ+ and 
Auj = A uj + . Therefore Gabor suggested to replace inequality (69) by the following 
modification: 


where 


Auj + At + > 1/2, (81) 

/ oo poo 

(t-i + ) 2 \f + (t)\ 2 dt, i + = t\u(t)\ 2 dt. 

-oo J — oo 

Of course, inequality (81) is correct. But now we meet another problem: what is 
the physical meaning of the quantities At + and t + , and how these new quantities are 
related to At and 1? It was shown by E. Wolf 54 that 1 + = t and A t + = At under 
the condition 


m = -}= [°° f(t)dt = o. 

V J — oo 


In such a case 


Aw+A t > 1/2, 


(82) 


(83) 


The proof of (83) was given by Kay and Silverman 55 without using the concept of 
analytical signal. We reproduce it (in a slightly modified form) in the case of t = 0, 
which can be always achieved by a shift of the origin of the time axis. 

Remembering that the Fourier transform of function tf(t) is idF/du> = iF'{u>), 
we can write (using the Parseval identity) 


pOO pOO 

(At) 2 (Auj + ) 2 = 4 / \F'(uj)\ 2 dui / (oj — w+) 2 \F(oj)\ 2 duj. 

Jo Jo 


(84) 


Using the Schwartz inequality (24) with /(w) = F'(co) and g(oj) = (oj — w + ) F(uj) we 
are led to the following chain of relations: 

r°° r°° d 

Auj + At > / (uj-uj+)[F'F* + (F')*F]duj = / (w - w+) — \F{uj)\ 2 duj 

Jo Jo 

poo 

{{uj-dJ + )\F(uj)\ 2 }\™ - \F{u:)\ 2 duj = 11/2 — |F(0)| 2 uJ+1 . (85) 

Jo 


However, not only the equality in ( |85| ) cannot be achieved (as was shown in [55]), but 
it gives in many cases the lower bound on the product Aco + At, which is significantly 
smaller than really admissible values. 

For example, let us consider the Gaussian signal with 1 = 0 [that does not satisfy 
condition (82)]: 

f(t) = (ct 2 /71-) 1/4 exp (-cr 2 ! 2 / 2 ) . 


F(uo) = ( 7 rer 2 ) exp (—w 2 /2ct 2 ) . 


( 86 ) 

(87) 


Then 


At = 


UO+ = 


Aw + = a 


7T — 2 


1/2 


27T 


a\/2 V 71 

so that for any a we obtain the same value of the product 

(AtAu+hauss = [(tt - 2)/(4tt)] 1/2 » 0.301 
whereas ill — 2|F(0)| 2 ciJ + | = i |1 — 2/7t| ss 0.18. 


( 88 ) 
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The exact lower bound for the product AtAw+ was found by Mayer and 
Leontovich 52 , whose results were reproduced in the book 56 . However, these 


studies were unknown in the West for a long time. Hilberg and Rothe 57 gave 
in 1971 a solution to the problem, very similar to the Mayer-Leontovich approach 
(except for the final stage). Then Borchi and Pelosi 58 gave a derivation, using a 
quantum-mechanical analogy with a double harmonic oscillator. We reproduce below 
the Mayer-Leontovich approach, following 56 . The reason is that this nice solution 


is still not well known. Besides, we give some details absent in 52,56 and compare 


their results with 57 58 


We assume that t = 0, since this equality can be achieved by a simple shift in 
time. It was proven in 52 that the minimum of product Aw+Af can be achieved 
under the condition t = 0 for even functions of time /(f). Then the spectral function 
F(u>) is also even and real (so that condition (82) is not satisfied). Abandoning for a 
while the normalization condition on function F(u), one can rewrite the right-hand 
side of (84) as a functional 0{F} of the function F(ui ): 

Q{F} = AN 2 /N 2 - AN^/Nq, 

where 


(89) 


A = 2 


[F'{u])} 2 duj , N m = 2 / u> m F 2 (to)duj. 


The variation 5Q must go to zero for the extremal function. This means that 
N$5Cl = N 0 (N 0 N 2 - Nf) 5A + N 0 A ( N 0 SN 2 - 2Ai<57Vi) 

+ A (mf - 2 N 0 N 2 ) 5N 0 = 0. 


(90) 


Let us introduce the notation a = A, b = N 2 and c = N\ for the extremal function. 
Imposing now the normalization condition N 0 = 1 (then c = w + ), we can represent 
equation (90) as 

6 j |a (w 2 — 2cw + 3c 2 — 26) F 2 (w)+(6 — c 2 ) [F'(w)] 2 |dw = 0.(91) 

This is a standard variation problem. The corresponding Euler-Lagrange equation 
reads 

(b - c 2 ) F"(u>) - a (w 2 - 2cw + 3c 2 -2 b) F = 0. (92) 

Multiplying equation ( [92] ) by F'(oj) and integrating from 0 to oo we get 

a (oj 2 — 2cw + 3c 2 — 26) -F 2 (cn)|^° = (6 - c 2 ) [E'(w)] 2 |^° =0, 

because F'(0) = 0 due to the parity of function F(u). Consequently 6 = 3c 2 /2. 
Introducing the notation 


Strain — [(AfAcC_|_) 


2 _ 2 
min\ — 


we obtain 




2 = - a b = - 


[F'{uj)] 2 duj / oj 2 F 2 (uj)dijj. 


3 3.. 

Then the Schwartz inequality yields the lower bound 

1 


2 

^75 


ujF' (ui)F(ui)du! 


v/12 


0.289. 


(93) 


(94) 
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In view of (88) we can conclude that 0.289 < p < 0.301 (the result p, « 0.3 was 
suggested also in [55]). To find the exact value of p, it is necessary to solve equation 


(|92|, which can be reduced to the known Weber equation [59 
d 2 F/dz 2 + (v + 1/2 - 2 2 / 4 ) F = 0 
with (taking into account the relation b = 3c 2 /2) 

v + 1/2 = \J ab/3 = p, z(ui) = 2 yfp(u>/c — 1). 


(95) 


Solutions to equation (95) that are limited for z —> oo are the parabolic cylinder 


functions D v {z). Consequently, the extremal frequency spectrum is given by the 
formula 

= V 1/2 (2 y/ifa/c - 1]) ■ (96) 

(Hilberg and Rothe [57 did not want to analyze this solution, choosing a longer way.) 
The extremal function F fl (uj) must satisfy the additional restriction F' (0) = 0. Since 
z(0) = —2 yfjl, the value of p can be found as the solution of equation 

K- r/a (~2 Vm) = 0, D' v {z) = dD v {z)/dz. (97) 

Using the integral representation of the parabolic cylinder functions with negative 
index v < 0 59 (in our case p ~ 0.3, so that v : 


- 0 . 2 ) 


one can 


t r 00 

D u(z) = yo-v / exp (-z 2 / 4 - zt - t 2 / 2) t~ v ~ 

r (-^) Jo 

reduce equation ([97]) to [52] 

/•OO 

/ exp (2pt — t 2 / 2) t~ p _1 / 2 (t — p)dt = 0, 

Jo 


l dt 


(98) 


where p = ^fp, . However, Mayer and Leontovich did not give a numerical solution to 
this equation, having stopped at this point. Our numerical solution of equation (98) 
gives the value 

/i = (Aw+Af) min = 0.29505306... (99) 

The value given by Hilberg and Rothe 57 was (one should divide their result 
1.180212... by the factor 4 due to the different definition of ‘uncertainties’) 0.295053, 
and the same value was given by Borchi and Pelosi 58 , who solved equation (971 
using another way. 


very close to the minimal possible value (99) (this was noticed in 56 


We see that Gaussian signals (|86|-( [87|) p ossess the value of product Acu + At that is 

Therefore it is 


interesting to compare Gaussian packets with the ‘minimal time-frequency uncertainty 
packets’ (96). We notice that the extremal frequency spectrum is more ‘flat’ near the 
point ijj = 0, since the consequence of equation (92) with b = 3c 2 /2 is the equality 
F"( 0) = 0. On the other hand, the extremal function decays slightly faster than 
the Gaussian signal for w —>• oo, since asymptotically D„(z) ~ z v exp (—z 2 /4) [59]. 
The comparison of the spectral functions (87) and (961 with identical values of is 


shown in figure |T[ Both the functions are normalized by their values at u> = 0, and the 
dimensionless variable w means in fact the ratio u/u + . Formally, this means that we 
take c = 1 in formula (96) and a = in (87). Figure [2] shows the time-dependent 
real Gaussian signal (86) and the real Fourier transform f^ft) of spectral function 
(|96[ as functions of the dimensionless variable Uj + t. [For u> < 0 we use the relation 


F(—u>) = F(uj).} Again, both functions are normalized by their values at t = 0. We 
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CO 

Figure 1. Gaussian spectral function F g (uj) J87j) and the ‘minimal uncertainty 
product’ spectral function ( |96| with equal dimensionless average values 

= 1. 



Figure 2. Time-dependent signals corresponding to spectral functions ( |87[ ) and 
( |96[ ) versus dimensionless time variable HJ+t with UJ .|_ = 1. 


see that f^{t) can assume negative (although relatively small) values. Similar figures 
were given in 


57 


but the ‘extra-flatness’ of function F^(uo) at its top was not noticed 













Energy-time and frequency-time uncertainty relations: exact inequalities 


17 


there. 

Approximate ‘minimal time-frequency uncertainty packets’ in the 
truncated Gaussian packets F(w) = exp [—(|w| — m) 2 /2cr 2 l were studied in 


55 


form of 
where 


0.3 for certain choices of parameters m and a. Spectral 

the 

points ±wi were considered in 


it was shown that Aw + Af 
functions in the form of a combination of two Gaussian functions centered at 

The minimal value of the product Af Aw+ 


60 


for 

such functions is 0.29535, but the functions f(t) have the form of Gaussian packets 
multiplied by cos(wif), so these functions also assume small negative values in some 
time intervals. Parabolic cylinder functions arise in the problem of minimization of the 
product At Aw in the case of band-limited spectral functions satisfying the condition 
F(w) = 0 for |w| > 61 . Some results of recent studies on the time-frequency 


problem can be found in 62 ■ 66 


6.1. Consequences for the AEAt problem 


Now we can return to the AEAt problem. Comparing formulas (70)-(741 with (10), 
(16), (59) and (60), one can see that the quantities e (60) and Ait (62) are equivalent 
to At (70) and Aw + ( f77| ) , respect ively{f[] Therefore, the lower bound of the product 
eAi t cannot be less than 0.29505A This boundary could be attained by replacing 
function [/(f)] 2 with [P(E)] 2 and [F(w)] 2 with 

Q*{t) = [x*(*)] 2 = G{D„_ 1/2 (2 rfi[t/c - l])} 2 , 


where the constant factor G ensures the fulfillment of the condition Q(0) = 1. But 
here we meet a difficulty: the energy distribution function P{E) must be non-negative, 
whereas the extremal function /(t) in the frequency-time problem can assume negative 
values, according to figure[2j Consequently, function Q*(t) cannot be a true non-decay 
probability. Therefore the minimal possible value of the producr eA\t is between the 
values 0.295fi and 0, 301 h (attained for the Gaussian non-decay probability function). 
The difference is rather small, but the question on the exact lower bound is still open, 
although, perhaps, it is not very important for applications. 


7. Mixed quantum states: the Eberly—Singh approach and its 
generalizations 



the framework of the quantum Liouville equation for the statistical operator (density 
matrix) p 

ihdp/dt = [H, p\. (100) 

Remember that the notion of statistical operator was introduced by Landau |69 and 
a few months later in a more elaborated form by von Neumann 70 . It is worth noting 


that equation (1001 is written in the Schrodinger picture. In the Heisenberg picture it 


reads dp/dt = 0, meaning that the statistical operator is the integral of motion (this 
fact was emphasized in 171 ). 

|| It is funny that the time variable in the frequency-time problem should be substituted by the 
energy in the time-energy problem, while the frequency should be substituted by time. 
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Following 18 let us introduce “the stationarity time” To according to the formula 

\2 


Tq" 2 = Tr 


(dp/dty 


( 101 ) 


Taking A = H and B = dp/dt = p in inequality (|4j), we have (assuming h = 1 in order 
to simplify formulas; the place of h can be easily recovered from the dimensionality 
analysis) 

(AH) 2 Tr (p/?) > 1 |w (p[H,p\) | 2 = \ |Tlr (p[H, [.H,p ]]) f . (102) 

Using the identity Tr B]C^j = Tr ([C, A]B ^ , we can rewrite Jl02|) as 


/ a2\ 1 


IV 

4^1 

[ Tr (' , )i 


(103) 

One can easily check the following identities: 

Tr (p^ 2 ) = Tr(p[p, H][H, p\) = Tr {p 3 H 2 - p 2 HpH) 

= * Tr (p 2 [ff, [H, p]]) = ±Tr (p 2 i) . (104) 

Let {|n}} be a basis where operator p is diagonal: 

P = ^2pn\n)(n\, 0<p n <l, ^2,p„ = 1. 

n 

A ^ a 2 

Then Tr [p 2 A) = YlnPn(n\A\n). Since the quantity Tr(pp ) is nonnegative for any 

statistical operator p, all mean values {n\A\n) are nonnegative, as well. Then 


Tr {p 2 A) < '^2,p n {n\A\n) = Tr (pA). 

n 

Consequently, relations ( |104| )-( [T05 ) lead to the inequality 

Tr (pf ) < ±Tr (p[H , [H,p]]) = ±Tr (j?) . 

Putting it into the left-hand side of inequality (103), we arrive at the relation 
(recovering the Planck constant) 


(105) 


(106) 


{ AH) 2 > h -Tr^). 


(107) 


Defining the energy “uncertainty” as A E = AH , we obtain the inequality 

AET 0 > h/y/2. (108) 

We can introduce a set of “stationarity times” T n according to the formula 


TU 2 = Tr 


p n {dp/dty 


They satisfy inequalities 

T n > T n _ 1 > • • • > T\ > V2T 0 , 


(109) 


( 110 ) 


following from (1051 and (106). In the case of pure quantum states, when p n = p, all 
these inequalities turn into strict equalities. Relations (110) and (108) result in the 
inequalities 


AET n > h, n > 1 


( 111 ) 
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(namely the special case of n = 1 was considered for the first time by Eberly and 


Singh 67 


It can be interesting to see explicitly, why the times To and T\ differ by the factor 

\/2 in the case of pure quantum states p = \ip)(ip\. The Schrodinger equation (with 

-2 

H = 1) \if>) = — iH\if) leads to the following operator p : 


= \iP){if\H 2 \-if){if\ - H\if){^\H\ip)(il>\ + H\tP)(m - \i/;)(t/,\H\il,)(4,\H, 


( 112 ) 


so that Tr (p ) = 2(A E) 2 . On the other hand, multiplying both sides of identity 

(112) by the operator |'i/’)(V’l from the left, one can see that the last two terms in (112) 
cancel each other. Then 

pp = \tp){if\H 2 \ip){if\ - 


and Tr {pf?^ = (AT) 5 


The advantage of inequalities (108) and (111) over the Mandelstam-Tamm 
inequality ([5]) is that the “stationarity times” T n are determined completely by the 
density matrix of the system alone, whereas the ‘duration’ A tA in ([5]) depends on the 
choice of an auxiliary operator A. A generalization of the Eberly-Singh inequality in 
the form 


T^ k AH > h/ 2 , 


rpLJK _ 

1 S ~ 


s = dp/dt 


(113) 


\(ds/dt}\ ’ 

was obtained in |72j; it was shown that Tg K < Ti, with an equality attained for pure 
quantum states. 


7.1. Influence of quantum purity 


Note that relations (108) and (111) become identities for any pure quantum state. 
Therefore it is interesting to know, in particular, how the lower bound of the 
product ATTo can depend on the quantum purity p = Tr(p 2 ). This problem was 
not considered earlier. A simple solution can be found in the special case of a 
harmonic oscillator and the initial Gaussian states (which remain Gaussian during 
the evolution). In this case it is convenient to use the description of quantum states in 

related to the complex matrix elements ( x\p\y) 


terms of the real Wigner function 73 


of statistical operator p in the coordinate representation as follows, 
W(q,p) = J dve lpv/h {q - v/2| p \q + v/2 ), 


W(q,p)dqdp/(2TrH) = 1. 


(114) 


(115) 


The advantage of using the Wigner function for (closed) quantum systems with 
quadratic Hamiltonians is that the evolution equation coincides with the classical 
Liouville equat ion, containing partial derivatives of the first order only. In particular, 
equation (100) with Hamiltonian H =p 2 /(2m) +moj 2 x 2 /2 can be transformed to the 
form 


dW 

~dt 


= tow q- 


dW 

dp 


p dw 

to dq 


(116) 
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Using the formula 74 77 

Tr (Ab) = J W A (q,p)W B {q,p)dqdp/(2Trh), (117) 

where W AB (q,p) are the Wigner-Weyl symbols of operators A and B, defined in the 


same way as in equation (114), we can rewrite definition (101) as 


T 0 -2 = J {dW/dt) 2 dqdp/(2Trh). 


(118) 


Let us consider the initial quantum state described by the most general Gaussian 

(119) 


Wigner function 78 


w{q ' p) = 7B exv V^D L 


0V> (q - qf +Vq(p- Pf - 2 a qp ( q — q) {p — p) 


where 


D = <T q a p - a 2 p = [(h/{2p)} 2 > H 2 / 4. 


( 120 ) 


The last inequality in (120) is the Robertson-Schrodinger uncertainty relation, which 
must hold for any pure 10 or mixed 78 quantum state. 

Calculating Gaussian integrals one can obtain the following formulas: 


rji —Z _ 

a n — 


H(B- 2Deo 2 ) 

8 T > 3 / 2 ; 


(at ) 2 = |-J(M 2 , 


2(AET 0 /H) 2 = (^Pj 


3/2 


1 +0J' 


,2 D-h- 


72\ 


B - 2 Doj 2 J 


( 121 ) 

( 122 ) 


B = m 2 ui 4 (a 2 + 2a q q 2 ) + (a 2 + 2a P p 2 ) /m 2 + 2w 2 ( a 2 p + 2a, 


qpq 


For pure Gaussian states the product ( fl22| ) equals unity for any choice of p aram eters 
(satisfying the restriction D = H 2 /4 under which the Wigner function (119) has 
sense 78j 79]). This shows a crucial difference between quantum pure and mixed 
states. For pure states, the stationarity time To is determined completely by the 
initial energy dispersion A E. On the contrary, the value of T 0 for mixed states with 
a fixed value A E 0 can vary in wide limits, including the value T 0 = oo for initial 
thermal states with muj 2 a q = a p /m = oj\Td > tno/2 (for w > 0) and q = p = a qp = 0. 
The minimal value of the product (122) for D > h 2 /4 equals (4 D/h 2 )^^ 2 . It is achieved 


if B —> oo. This can happen, for example, for highly squeezed thermal states (a q —l 0 
or a p —> 0 with D = const, so that a p —> oo or a q —>• oo, respectively) if q = p = 0, 
strongly correlated mixed states ( a qp —t oo with D = const and q = p = 0) or for 
strongly shifted states (q —> oo or p —>■ oo). Thus we arrive at the inequality 

A£T 0 > H(2p 3 y 1/2 (123) 

which holds for any Gaussian state. The equality sign takes place for free Gaussian 
packets with u> = 0. We see that pure quantum states are ‘more fragile’ than mixed 
states with the same energy dispersion: j , ( pure ^ < Tq 


Inequalities (111) can be strengthened, if one takes into account the relation 


Tr (AB) < Tr(^4)Tr(R), which holds for Hermitian positively definite operators. 
Introducing the ‘higher order purity’ p n = Tr(p"), we obtain the inequalities 

A ET n > AETop- 1 ' 2 > H(2p n )~ 1 / 2 , 


(124) 
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A ET n > AETuj-l{ 2 > n^[ 2 , n> 2. 


(125) 


It would be interesting to generalize inequality (1231 for arbitrary (including non- 
Gaussian) mixed quantum states. Probably, a simple function /i -3 / 2 can be replaced 
by some more complicated dependence on the quantum purity. Such a hypothesis 
follows from a similar example related to the Heisenberg-Kennard-Weyl inequality 
([3]). The product AxAp equals h/(2p) for thermal states of a quantum oscillator, but 


the general ‘purity-dependent uncertainty relation’, obtained in 19 80 , has the form 


AxAp > 2, where function $(/x) has rather complicated explicit expression, 

going asymptotically to <l>(/z) « 8/(9/x) for p <C 1. Perhaps, simple inequalities (124) 


and (125) can be also strengthened and generalized. 


7.2. Remarks on open (dissipative) quantum systems 

It is worth noting that both the Mandelstam-Tamm inequality © and inequalities 
(123) (1251 are essentially based on the fact that the evolution of a quantum system is 


unitary , being described by the linear Schrodinger equation or its consequence (100), 


as well as on identifying the energy with the (Hermitian) Hamiltonian. For other kinds 
of evolution equations, e.g., for the ‘standard master equation’ [J] 


dp 

dt 


N 

£ 7 , 

i=i 


2AjpA\ — A\Ajp - pA\A 3 


- 


(126) 


describing open (dissipative) quantum systems, inequalities (|5| and (123)- (125) do 
not hold. This can be easily seen already in the simplest example of the relaxation of 
a quantum harmonic oscillator due to the contact with a thermal bath at zero absolute 
temperature. Suppose for simplicity that the initial density matrix was diagonal in 
the Fock basis. Then it is known that all off-diagonal elements will remain zero with 
the course of time, whereas the evolution of the diagonal elements (level population 
probabilities) p n = p nn = (n\p\n) is given by the equation [following from ( |126[ ) in the 
case of j = 1 and A\ = a - the boson annihilation operator] p~| 

dp n /dt = 27 [(n + l)P„+i - np n ] , n = 0 , 1 ,... (127) 

For the initilal Fock state | M) with M ^ 0, the initial energy dispersion A E = 0, 
while 


T 0 2 = £(rfPnM) 2 = (2 7 M) 2 

n—0 


at the initial instant t = 0. Consequently, AETq = 0 in this case. 


It was suggested in 87 that the irreversible evolution can be compatible with the 


time-energy uncertainty relation for some class of nonlinear evolution equations for 
the statistical operator p, with a suitable definition of the characteristic evolution 
time. However, the situation is different for many nonlinear generalizations of 
the Schrodinger equation, which were proposed in attempts to describe effects of 

K This equation was derived by several authors in 1960s |8l|)83| ; but nowadays it is called frequently 
as ‘Lindblad equation’ after paper 84]. 

+ This equation was derived for the first time by Landau 69] and generalized later by several authors 
in 1960s: see, e.g., [85||86]. 
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dissipation in quantum mechanics in terms of the wave function (see, e.g., review 
The concrete example of equation proposed by Kostin l8fT 


QJj A 

= H e fijj = V 2 ?/; + V(x)ip - ih'y 


In I t I - 


In it 


if 


(128) 


y if*J \ \if * 

was studied in 90 for the harmonic oscillator potential in one dimension V(x) = 


muj 2 x 2 /2. Equation (128) admits exact solutions in the form of time-dependent 


Gaussian wave packets. Analyzing such solutions and calculating the time-dependent 
characteristic time To(t) pit , it was discovered that equality 


(AH ef ) 2 T 2 (t)=h 2 /2 


(129) 


is fulfilled for any solution at any time instant. However, the physical meaning of 
the variance (A H e f) 2 is not clear, since the ‘effective Hamiltonian’ H e f depends on 
the wave function. On the other hand, if one uses the natural definition of ‘energy 
uncertainty’ A E as the dispersion of the standard energy operator + K(x) 

calculated with respect to the time-dependent solution if(x,t), then AE(t)To(t) ss 
hco/(2\/2'-f) in the over-damped case |y| w, so that this product can be made as 
small as desired. 


8. Speed of quantum evolution 

The energy-time uncertainty relations are closely connected with the problem of 
finding the minimal time to achieve some reference subspace in the Hilbert space, 
starting from a given initial state; in particular, to reach some state orthogonal to 
the initial one. This subject was studied by many authors, so we give here the 
most relevant results and references only. The first papers in this direction seem 
In particular, considering the unitary evolution of pure quantum states 


to be 91-96 


governed by time-dependent Hamiltonians H(t ), the following inequality was obtained 


91 92 


(AE)At± > ttH/2. 


(130) 


Here At± is the minimal time interval necessary to achieve a state orthogonal to 
the initial one, whereas (A E) is the time-averaged uncertainty in energy during this 


interval. The discussions of (130) and further generalizations (in particular, for mixed 
quantum states) can be found in 28 97 105 . We bring here only one of results 
obtained in 196], generalizing ^ and (|9|: 

sin (6 - h t ) < \{cp\ip t )\ < sin (5 + h t ). (131) 

Here | iff) is the time-dependent pure quantum state originating from l^o) at t = 0, 
| (p) is some fixed quantum state, <5 = arcsin|(</j|'i/>o)| and 

ht=h ~ 1 [ min {A Lp H(s), A^ 0 H(s)} ds, 


where A V H = 


(v\H 2 \y) - (vWWY 


If \tp) = |^o), then 


1/2 


|(V ; o|V’t)l > cos [H / A ^ 0 H(s)ds 


(132) 
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If | ip) _L IV'o), then 

l(¥#t)| < sin (h t ). 


(133) 


This means, according to [96], that a state orthogonal to the initial one cannot be 
populated too fast, and that a fast transition to an orthogonal state can happen if 
only both energy uncertainties, A ip H(s) and A ^ 0 H(s), are large enough. 

8.1. Margolus-Levitin inequalities and their generalizations 

In the inequalities considered in the preceding sections the quantity A E was related 
somehow to the energy dispersion , defined by various ways. A new kind of relations 
containing the difference between the mean energy and the minimal possible value of 


energy for the given physical system was discovered by Margolus and Levitin 106 


Following them we assume in this subsection that E = 0 is the minimal possible value 
of energy. Then using the inequality 

cos(a;) > 1 — (2/7r)[a: + sin(a;)] (134) 

one can write the following relations for the real and imaginary parts of nondecay 
amplitude x(t) ([T6[: 


Re[x(t)] = [ cos (Et/h)P(E)dE > 1 — —^~t + — Im[x(t)]- (135) 

Jo nn tt 

Let Tj_ be the time it takes for |^>(0)) to evolve into an orthogonal state, so that 
XCTO = 0. Then Re[x(Tj_)] = Im[x(Tj_)] = 0. Consequently 

{E)T± > nh/2. (136) 


Applications of the quantum speed limit inequalities like (136) to the areas of 


quantum computation, quantum information and quantum metrology were discussed 

in 107 111. 

For example, Luo and 


Generalizations of (136) were obtained in 112 -120 


Zhang 114 introduced the quantity 

T a = infjf > 0 : Q(t) = a} , 
and proved that (here H = 1 for simplicity) 

m ^ _ ( 1 - \A*(! + VA 2 )' 

a - ^ y 2 (EP) 

for 0 < p < (Tr/2)y/(lJa)--T. In particular, 

To = T± > tt(2 {E p ))~ 1/p 


0 < a < 1 


i /p 


(137) 


(138) 


(139) 


(the quantity To used here is different from that defined by equation (101) ). The 


main inequality used for the derivation of (1381 is 


cos(a;) + (2p/tt) sin(.'r) > 1 — 2(x/tt) p . (140) 

For p = 1 inequality (138) can be easily derived from ( 135| ), if one takes into 
account that for the fixed value Q = \x\ 2 = a one can write Re(x) = v / ' acos ( < / > ) 
and Im(x) = i/^ s i n (^ ) ) with some phase <j>. Since the maximal value of function 
/((/>) = cos {(j>) — asin(</>) is \J 1 + a 2 , one has 

1 - < Re[x(i)] - —Im[x(t)] < y/a (i + 4 /tt 2 ). 

7 Tti 7r 
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Consequently 


2 (E)T a >nh(l- y/a{ 1+4/tt 2 )) . 


The following chain of generalized inequalities connecting the time To with higher- 
order moments of the energy distribution function was derived in 117 : 


{(E ~ {E)) 2n )T 2 
(2 n)\H 2n 


> 


n 1 ziy 

£ — 

k =0 


— k 


{{E-{E)) 2k )T 2 


(2 k)\h 2k 


Note however that the first inequality of this chain (for n = 1) T 0 AE > h\J 2 is weaker 
than the M-T bound (130), since \[2 ~ 1.41 < 7r/2 ps 1.57. The states that optimize 
inequality ( |136 ) or its generalizations (related to the problem of time-optimal quantum 
evolution) were considered in 121 -127 . 

The situation is more complicated for mixed quantum states, especially in the case 
of non-unitary evolution (for open quantum systems, for example) [127 134 . First of 
all, one has to generalize simple formula (10), which seems a natural definition of the 
‘non-decay probability’ for pure states. The authors of papers [128 131 considered 
the ‘relative purity’ 

Q 1 (t)=Tr[p(t)p(0)]/Tr[p 2 (0)], (141) 

whereas the authors of fl!2l[Il3l[IT8l[l29l[l30l used the so called Bures-Uhlmann 
fidelity 135 137 


Q-2(t) = 


Tr 


1/2 - - 1/2 
Po PtPo 


(142) 


where operator p 1 / 2 is defined as the unique positively sem idefinite Hermitian operator 
satisfying the relation (p 1 / 2 ) = p. One of results of 


131 


was the inequality 


T* 


(1 - Q*)Tr (p; 


/Tr 


(&Po) 


(143) 


where C is the generator of evolution defined according to the equation dp/dt = Cp and 
T* is the time necessary to diminish the value of Q\ from unity to Q*. Obviously, the 
initial energy dispersion is not the decisive factor in the case of non-unitary evolution. 
Qualitatively similar results in paper [130 were expressed in less explicit form in terms 
of the quantum Fisher information F(t) = Tr p(t)L 2 (t ) , where Hermitian operatot 

7(f) is defined implicitly through the evolution equation dp/dt = ('Lp + pL^j /2. 
Counterexamples showing that there is no quantum speed limit for non-Hermitian 
Hamiltonians or a nonunitary evolution were given in [129,138 . Other inequalities 
relating the ‘distinguishability time’ with the quantum Fisher information of p on the 
unitary path generated by H (and with AH for pure states) were obtained in 139 . 


8.2. Arrival and tunneling times 

Another large field of research connected with ETUR is the problem of arrival time 


15 140 146 . For example, it was shown in 141 that the arrival time of a free 


particle to some fixed point in space cannot be deterined with the precision better 
than At arr i v > h/(E), where (E) is the mean initial kinetic energy. Note that this 
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145 


relation resembles (1361. This result was generalized to the case of arbitrary potentials 


be found in 22 


A comprehensive review of studies on the arrival time problem until 2000 can 
while for more recent publications one can consult e.g. 


The concept of ‘quantum delay time’ (in scattering) was analyzed in 160 162 


The problem of tunneling time (or traversal time) was considered in 161 163-177 


147 159 


Connections between the energy-time uncertainty relations and the quantum Zeno 
effect 178 were discussed in references 179]-183. On the other hand, the uncertainty 
relation in the non-strict form 0 was used in 184] to explain a possibility of anti-Zeno 
effect. 


9. Problem of time operator 

Inequalities 0 or 0 could be derived immediately from the commutation relation 
[H,T\ = -ih, (144) 

if the time operator T existed. However, it was noticed by Pauli as far back as in 


1926 that no Hermitian unbounded operator satisfying (144) can exist for an arbitrary 
Hamiltonian H 


185 


This is connected with the specific property of energy spectrum 
of physical systems: since the spectrum of the time operator must be undoubtedly 
continuous and unbounded, the same properties must also be inherent in the spectrum 
of the Hamiltonian. However, energy spectra of the majority of physical systems are 
bounded from below; in addition, they can be discrete. Nonetheless, many people 
tried to find some surrogates of the time operator for various specific systems or in 
some restricted sense. For example, such a problem arises in scattering theory, if one 
tries to define the duration of a collision 


15 186 


In the simplest case of a free particle Hamiltonian H 0 = p 2 /(2m), remembering 
the classical formula x = pt/m + const , one can define a formal ‘time operator’ 
as fl4l 187 188 


T 0 =m(p 1 x + xp *) . (145) 

In the three-dimensional case Razavy |l89 considered the ‘time operator’ in the form 


Tr = y {fp r +p r f)H r \ 

where p r and H r are the radial momentum and radial part of Hamiltonian, 
respectively. 

Although the operators Hq and T 0 formally satisfy relation (144), it is clear that 
operator (145) has a lot of defects. First of all, it contains a singular operator p _1 . 


In addition, although operator (145) in the energy representation can be reduced to 


the form To = ihd/dE , nonetheless it is non-Hermitian in the space of functions used 
in physics usually. Indeed, the hermiticity condition demands the eigenstates ip{E) to 
form a complete set in the semiaxis E > 0, to turn into zero at E = 0 and to be closed 
with respect to the operation d/dE. Such a set of functions does not exist, as was 
pointed out in 15]. Therefore when dealing with operators like (145) one should either 


ignore their unpleasant properties or define the class of admissible wave functions in 
a special manner (see also [190,191 in this connection). The proposal to abandon the 
idea of unbounded time and to use instead some bounded time operators satisfying 
(11441) was developed in 1921. 


Note that if the question about the existence of ‘time operator’ is understood 
in a restricted sense, e.g. as the question of finding an operator satisfying relation 
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(144) for a given Hamiltonian, then such an operator can be found probably for any 
In particular, the prescription for how to do this for one-dimensional 
But the real problem consists in finding conditions 


193 


Hamiltonian, 
systems was given in 
under which the formal ‘time operator’ proves to be Hermitian (more precisely, self¬ 
conjugate). Even more important is the possibility of physical interpretation of such 
an operator as true time operator. For example, the problem of the energy spectrum 
boundedness does not exist for the Hamiltonian Hp = p 2 /(2m) — Fx describing a 
particle moving in a uniform potential field. The spectrum of Up is continuous 
and extends from — oo to o o. The operator Tp = p/F 189 190 194 is Hermitian 
and satisfies equation (144), but what is its relationship to time? Moreover, Tp 
does not commute with the coordinate operator x. Thus, if it were the real time 
operator, this would mean that it is impossible to determine simultaneously the 
coordinate of a particle and the time at which this particle passes through the point 
with this coordinate. For these reasons the difficulties with physical interpretation 
force us to treat most ‘time operators’ constructed thus far as purely mathematical 
artificial constructions without true physical meaning. The worst feature of such ‘time 
operators’ is that they are not universal. Instead, they should be adjusted each time 
to the concrete Hamiltonian: compare operators T 0 and Tp. 

Nonetheless attempts to construct operators resembling time in some restricted 
sense continue. One direction is to extend the Hilbert space, thus removing the 
problem of lower bound of the effective Hamiltonian 195 196 . A time-like operator 
for the harmonic oscillator was constructed in 197 and for a singular oscillator 
More general cases of systems with discrete energy spectra were considered 
In this connection, the entropic time-energy uncertainty relation was 


m 

in 


198 


199 


201 


introduced in 201 


One may suppose that the time operator could arise naturally in the relativistic 
quantum mechanics: if operators exist for p = 1,2,3, then operator Xq must 
exist as well due to the relativistic invariance. Different approaches to constructing a 
relativistic time operator can be found, for example, in studies 202 ■ 205 . However, no 
unambiguous and generally accepted results were obtained in this field. Perhaps the 
reason for this failure is that in the relativistic case not only does the time operator 
not exist, but well-defined operators of the coordinates do not exist either. Thus, the 
time operator can apparently be introduced with the same degree of conventionality as 
the coordinate operator. For example, such a conventional operator was constructed 
by analogy with the known Newton-Wigner coordinate operator 
(in the momentum representation for a free particle): 

d Ed E 


206 in study 207 


T rP i — —ih [ 

i dE 


d 

-7T- + 


c 2 p dp 2 c 2 p 2 


(146) 


Then AEAT > H/2, but operator T re i does not commute with the momentum 
operator: 


Trel 7 P 

p E 

= in , 9 , 

■i J 

L 

c z p z 

J 


ihE 

c 2 p 


(147) 


The consequence of these relations is the rigorous version of the known Landau-Peierls 


inequality 208 


ATAp > *) — ^/(2c), 


(148) 
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which relates the accuracy of the measurement of momentum to the duration of the 
measurement. Other approaches were suggested or considered in 


209-216 


For other studies related to the problem of time operator or estimation of time 
displacements (durations) we can cite references 2l||217 - 258 where different solutions 
were proposed or criticized. 


10. Measurements of energy and time 

The most interesting (and the most controversial) field of applications of inequality 
|T]) seems to be quantum measurement theory. However, the existence of a relation 
like ([I]) in this field is still under question. Note that this question apparently did not 
exist for the majority of creators of quantum theory. For example, a number of famous 
thought experiments confirming relation 0 were proposed by Bohr in his discussions 


with Einstein during the Solvay congresses at the end of 1920s 259 . Since that time his 


examples have been widely reproduced in textbooks on quantum mechanics. However, 
Bohr’s arguments contained a weak spot, since he did not give a strict and clear 
definition for the quantities A E and At. Therefore all his relations are no more than 
estimates, so that inequality ([!]) is valid only to within an order of magnitude. For this 
reason authors of subsequent studies attached different meanings to the notion of the 
‘uncertainty in energy’. For example, this quantity was defined either as the accuracy 
of energy measurements, or the difference of energy values at two time moments, or 
the uncertainty of this difference, or the uncertainty of the ‘increase’ in the system’s 
energy during the time of measurement, and so on. Moreover, it was not quite clear 
what is understood by the term ‘energy’: whether it is kinetic energy (as soon as free 
particles only were discussed in the majority of papers devoted to this subject, the 
corresponding energies could be nothing but kinetic energies), or the total energy, i.e., 
the Hamilton function. In addition, no hint was given to a method of calculating the 
quantity A E for a given wave function or density matrix. The same can be said about 
the quantity At. Different authors defined it as ‘the duration of the measurement 
process’, or ‘the uncertainty of the time moment at which the change of energy of the 
system has taken place’, or something else. Needless to say, these quantities were not 
defined with mathematical accuracy. It is worth mentioning in this connection that, 
for example, such a quantity as ‘the duration of measurement process’ can be defined 
quite unambiguously in the special ideal case when the interaction between the system 
under study and the ‘measuring device’ can be turned on and off abruptly at some 
given instants of time. But the very existence of such abrupt interaction switching has 
been denied categorically by Fock in his papers (260 |261| devoted to the defense of 


relation ([!]). If, however, the interaction is switched on smoothly, then a constructive 
definition of its duration, absent in most old papers, is needed. 

For these reasons, the belief in the universal validity of relation 0 for all 
measurement processes has dwindled away, especially when some measurement 
schemes contradicting inequality ([!]) were proposed. The first investigations in this 
direction were performed by Aharonov and Bohrn 
other studies, 


e.g. 


262-266 


14 , which were followed by many 


For the sake of definiteness we put the following meaning to the quantity A E: it 
is an estimate of the dispersion of energy values before and after the measurement. In 
this case the question of the validity or invalidity of inequality ([Tj) is equivalent to the 
question of the validity of one of two statements: 

1. In any measurement continuing for a finite interval of time any quantum 
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system inevitably changes its state, so that the characteristics of the system before 
and after the measurement (for example, its energy) are inevitably different, and the 
dispersion of energy values before and after the measurement satisfies inequality 0 - 
where At is some characteristic duration of the measurement. Specifically, arbitrarily 
fast measurements of any characteristic of a quantum system without changing its state 
(more precisely, changing the state by some arbitrarily small amount) are impossible 
in principle. 

2. The possibility of arbitrary fast ‘nondemolition’ measurements does not 
contradict the principles of quantum mechanics (i.e., situations are possible when 
the product AEAt is as close to zero as one desires). 

Note that the problem consists in the reproducibility of the results of the 
measurement, because a single measurement of energy can obviously be performed (in 
principle) with an arbitrary accuracy (otherwise we would arrive at a contradiction 
with the postulates of quantum mechanics, since the energy observable corresponds 
to a Hermitian operator). The measurement, as a rule, changes the state of a system, 
but it is not clear whether such a situation takes place for any measurement, or 
whether somebody could invent such a scheme that the measurement would be exact 
and instant and the stage of the system would not be changed. Today this problem 
is important not only from an academic (or abstract) standpoint, but also from a 
practical (experimental) point of view, because it has an immediate relation to, for 
example, constructing detectors of gravitational waves which ought to work in the 
regime of the so-called quantum nondemolition measurements (QNDM) [268 -27 1]. 

In favor of the first statement we refer to classical papers (for example 208 259]) 
and the majority of textbooks on quantum mechanics. However, as has already been 
mentioned, only some examples have been discussed in that studies,, but no rigorous 

272]). On the other hand, it would seem that the 


proof has been given (see e.g. 
existence of a single counterexample is sufficient to refute inequality ([I]). Although 


several such counterexamples were given in studies 14 262 - 267 , the situation is 
nonetheless not so simple. There is a difference of principle between the examples 
in studies 208 259 and those in studies 14 262 - 267 . The authors of papers ‘in 


defense’ of inequality |l]) considered various measurement schemes which could be 
realized physically, if only in principle, in the real world. On the contrary, the papers 
like 14 262-267 [so to speak ‘against’ relation 0 ] deal with some mathematical 


models, i.e., their authors, as a rule, invent various expressions for the interaction 
Hamiltonians of the system under study with ‘measuring devices’ and investigate the 
properties of solutions to the corresponding equations. However, the way to realize 
these Hamiltonians in real physical experiments remains unclear. For example, a 
typical Hamiltonian of this type is as follows 14 263 266 


H = p 2 x /2m + p 2 v /2M + 7 (t)yp x 


(149) 


where a particle with mass to and coordinate x represents the system under study, 
while a particle with mass M to (in order that it could be considered quasi- 

classically) and coordinate y represents ‘the measurement device’. The interaction 
between the two systems is described by the term 'j(t)yp x , where 7 (f) is some known 
function of time disappearing for t —> ±oo. Fock [ 260J severely criticized Hamiltonians 


like (149) on the ground that they contain explicit functions of time. According to 
Fock, in this case we have to take into account the quantum nature of fields which 
are created by means of such time-dependent interactions (since, as was shown by 
Landau and Peierls 12 


only static fields can be considered, generally speaking, 
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classically, i.e., without taking their quantum nature into account). However, it 
was shown in 266 that we can take into account the quantum nature of the ‘field’ 


creating the interaction, believing 7 (f) to be the average value of some bosonic Held 
in a coherent state. Then the conclusion about the possibility of an arbitrarily fast 
and nondemolition measurement remains valid. (Strictly speaking, certain divergent 
integrals arise, as occurs frequently in the quantum field theory, due to the presence of 
an infinite number of field modes; but cutting off the frequency spectrum from above 
by the frequency u] max ~ (At) -1 , which seems quite reasonable and corresponds to 
the usual methods of quantum field theory, we obtain corrections tending to zero as 
At —> 0.) Thus we believe that it is more important to pay attention not so much to 
the function 7 (f) in formula (1491, but to the structure of the interaction Hamiltonian 
Hint ~ VPx■ how one could implement such an interaction in real experiments?! 

In other words, the answer to the question, which of the two statements 
formulated above is correct, depends in a decisive way on the answer to the following 
question: whether or not arbitrary Hermitian interaction Hamiltonians can be realized 
in our physical world? If the answer is ‘yes’, then statement 2 is correct; if the answer 
is ‘no’, then most likely statement 1 is valid. 

In summary, we will formulate our vision of the problem. It is based on 
distinguishing between quantum theory as a fundamental science studying the world 
around us (and not yet completed) and quantum mechanics understood as a 
phenomenological and nonrelativistic model which is (although logically closed and 
completed) suitable for describing a limited class of phenomena: see, for example, 
the discussion of this difference in 273 . With such an understanding of quantum 


mechanics as a mathematical scheme the second statement apparently seems correct. 
As to real cases, i.e., experiments which could be performed in the physical laboratory, 
inequalities like (|T]) must hold under the assumption (made implicitly by the founders 
of quantum mechanics) that the interaction Hamiltonians realizable in our world 
are those containing functions of the coordinates only (this was clearly formulated 
by Landau and Peierls [208] : they mentioned that if one could use any interaction 
Hamiltonian H ( p , x) , then the momentum of a particle could be measured for an 
arbitrarily short period of time without any change in its velocity). Still we will again 
note that, first, the exact values of the right- and left-hand sides of inequality ([l]) 
for such ‘real’ experiments have not been defined rigorously in many old papers, and 
secondly, the ‘derivation’ of this inequality occurs frequently inherently contradictory, 
especially in the cases when people try to prove 0 as a consequence of inequality ©• 


Such contradictions and ambiguities were analyzed in detail in 266 . For other studies 


on this subject we can recommend |8,9, 2lQp74p83] . Among many interesting results 
obtained in this area we bring two inequalities proposed by Vorontsov ; 265| (although 
they were not proven rigorously): 


AE + Ae > 


h 

At ’ 


AEAe > 


( 


h 

\2At 


(150) 


Here AE is the error of the measurement of the system’s energy during time interval 
At and Ae is an increase of the uncertainty of the energy of the apparatus (the 
disturbance of the energy of the apparatus). 
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We have shown that the area of energy-time uncertainty relations continues to attract 
attention of many researchers until now, and it remains alive almost 90 years after 
its birth. It received a ‘new breath’ in the past quarter of century due to the actual 
problems of quantum information theory and impressive progress of the experimental 
technique in quantum optics and atomic physics. It is impossible to describe various 
applications of the ETUR to numerous different physical phenomena in this mini¬ 


review. Instead, we conclude with a few relevant references 284-300 
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